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Abstract 

Aircraft  pitch  response  is  a  vital  element  of  piloted  vehicle  flying  qualities.  There  has 
been  controversy  over  both  the  form  and  the  substance  of  the  requirements  for  short-term  pitch 
response.  Currently,  MIL-STD-1797A  offers  six  different  methods  for  evaluating  short-term 
pitch  response.  These  six  methods  often  give  conflicting  results.  Two  methods  are  analyzed 
in  this  thesis  -  the  Neal-Smith  criteria  and  the  uspTe2,(sp,  t()  criteria.  Domains  from  both 
criteria  are  mapped  into  each  other  identifying  regions  of  conflict  and  regions  of  agreement. 
Parametric  studies  are  performed  and  evaluated  for  trends.  Further,  a  real-time  analysis  tool 
for  evaluating  these  methods  is  developed. 
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A  COMPARISON  OF  THE  NEAL-SMITH  AND  ujspTe2 ,  C.sp,  re  FLYING 

QUALITIES  CRITERIA 


I.  Introduction 

The  short-term  pitch  response  of  an  aircraft  is  a  vital  element  of  flying  qualities  (4:171). 
MIL-STD- 1 797 A,  Flying  Qualities  of  Piloted  Aircraft,  offers  six  different  methods  for  evaluat¬ 
ing  short-term  pitch  response.  There  are  strengths  and  weaknesses  for  each  method  depending 
on  aircraft  classification  and  flight  phase.  All  six  methods  have  been  maintained  because  the 
short-term  pitch  response  characteristics  are  universally  regarded  to  be  extremely  important- 
so  important  that  controversy  over  both  the  form  and  the  substance  of  requirements  still  exists 
(4:171).  MIL-STD- 1 797 A  provides  a  large  amount  of  guidance  for  determining  the  appropri¬ 
ate  method.  Still,  one  must  decide  which  of  the  six  methods  to  include  in  a  specification.  The 
simple  answer  of  including  all  six  leads  to  conflicting  results. 

Two  methods  are  analyzed  in  this  thesis  —  the  u>spTo2  Xsp,Te  criteria  and  the  Neal- 
Smith  Criteria.  The  uspTg2,Csp,Te  criteria  provides  a  relatively  simple  way  to  predict  flying 
qualities.  When  the  airplane  data  from  (1)  is  used,  the  criteria  has  a  confidence  level  of 
about  82  percent  (4:206).  Still,  the  uspTe2,(sp,  t0  criteria  only  deals  with  the  aircraft  as  an 
open-loop  system.  A  modified  version  of  the  Neal-Smith  criteria  uses  a  mathematical  model 
of  the  pilot  to  determine  the  pilot- vehicle  system  response.  This  closed-loop  criteria  is  known 
as  the  Pilot-in-the-Loop  criteria.  The  original  Neal-Smith  criteria  is  more  restrictive  than  the 
Pilot-in-the-Loop  criteria  in  that  pilot  phase  contribution  is  taken  into  consideration.  Other 
than  that,  the  two  criteria  are  essentially  equivalent.  Both  will  be  compared  to  the  open-loop 
-  Cspi  r#  criteria. 

One  method  for  comparing  the  iospTg2,(sp,  t6  and  Neal-Smith  criteria  is  to  map  one 
criteria  into  the  domain  of  the  other  criteria.  This  identifies  regions  of  conflict  and  regions 
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of  agreement.  Ultimately,  flight  testing  of  various  points  from  the  regions  of  conflict  could 
determine  whether  a  criteria  needs  to  be  modified.  The  answer  may  be  to  combine  the  criteria 
somehow  to  form  a  new  criteria.  The  motivation  for  this  type  of  research  comes  directly  from 
MIL-STD-1797A: 

The  large  amount  of  this  guidance  reflects  the  importance  of  short-term  pitch 
response,  the  high  attention  it  has  been  given,  and  the  great  need  for  further  study 
to  derive  a  clear-cut,  generally  applicable  set  of  requirements.  (4:172) 

1.1  Background 

Flying  qualities  is  that  discipline  in  aeronautical  engineering  that  is  concerned  with 
basic  aircraft  stability  and  pilot-in-the-loop  controllability.  ‘Flying  qualities,’  ‘stability  and 
control,’  and  ‘handling  qualities’  are  three  terms  which  are  generally  considered  synonymous 
(13).  To  prevent  confusion,  the  following  definitions  will  be  used.  Both  the  US  Air  Force  and 
Navy  agree  that  “flying  qualities  are  those  stability  and  control  characteristics  which  influence 
the  ease  of  safely  flying  an  aircraft  during  steady  and  maneuvering  flight  in  the  execution  of 
the  total  mission”  (10).  Edkins  defines  stability  as  “...the  tendency  or  lack  of  it,  of  an  airplane 
to  fly  with  wings  level”  and  control  as  “...steering  an  airplane  on  an  arbitrary  flight  path” 
(5).  Cooper  and  Harper  define  handling  qualities  as  “...those  qualities  or  characteristics  of 
an  aircraft  that  govern  the  ease  and  precision  with  which  a  pilot  is  able  to  perform  the  tasks 
required  in  support  of  an  aircraft  role”  (2).  Figure  1.1  shows  how  flying  qualities,  stability, 
control,  and  handling  qualities  are  related. 

Stability  and  control  analysis  deals  with  the  interaction  of  the  aerodynamic  controls 
with  the  external  forces  and  moments  on  the  aircraft  (3).  That  is,  stability  and  control  analysis 
primarily  deals  with  open-loop  systems  that  are  still  in  the  design  phase.  On  the  other  hand, 
handling  qualities  assessment  deals  with  the  pilot  and  aircraft  performing  as  a  closed-loop 
system.  This  assumes  the  aircraft  is  already  built.  The  six  methods  offered  by  MIL-STD- 
1797A  are  used  to  predict  handling  qualities  while  the  aircraft  is  still  in  the  design  phase. 
If  a  mathematical  model  is  used  for  the  pilot,  a  closed-loop  stability  and  control  analysis  is 
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Figure  1.1  Flying  Qualities  Breakdown 


possible  (e.g.  the  Neal-Smith  criteria).  Since  aircraft  perform  a  wide  variety  of  maneuvers 
and  vary  in  size,  some  type  of  grouping  is  necessary  before  analysis  can  begin. 

Aircraft  are  procured  to  perform  defined  missions.  The  class  designation  is  used  to 
help  particularize  the  requirements  according  to  broad  categories  of  intended  use  (4:76). 
The  intended  use  of  an  aircraft  must  be  known  before  required  configurations,  loadings  and 
operational  flight  envelopes  can  be  defined  (4:76).  MIL-STD-1797A  identifies  four  classes  of 
aircraft  as  outlined  in  Table  1.1. 


Table  1.1  Aircraft  Classification 


Class 

Description 

Example 

I 

small  light  aircraft 

T-41 

II 

medium  weight  aircraft  with 
low-to-medium  maneuverability 

C-21 

m 

large,  heavy  aircraft  with 
low-to-medium  maneuverability 

KC-10 

IV 

highly  maneuverable  aircraft 

F-16 

This  research  examines  two  different  aircraft.  The  first  is  Calspan  Corporation’s  variable 
stability  Lear  jet  24  (Fig.  1 .2).  The  second  is  the  Variable  Stability  In-flight  Simulator  (VISTA) 
which  uses  the  F-16D  as  its  host  aircraft  (Fig.  1.3).  The  Lear  jet  can  be  considered  a  Class 
II  aircraft.  The  F-16  VISTA  can  be  considered  a  Class  IV  aircraft.  Using  variable  stability 
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Figure  1.3  Picture  of  F-16  VISTA  (17) 
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aircraft  makes  it  possible  to  flight  test  the  analytic  results  of  this  thesis  in  the  future.  The 
cockpit  environment  can  be  changed  to  match  that  of  another  aircraft.  Since  an  actual  aircraft 
is  being  flown,  it  provides  a  degree  of  realism  not  present  in  ground  simulation.  As  the  pilot 
moves  the  controls,  he  or  she  experiences  the  real  flight  motions,  accelerations,  and  handling 
qualities  of  the  simulated  aircraft.  This  realism  gives  the  pilot  a  higher  level  of  confidence 
when  determining  a  handling  qualities  level. 

In  this  thesis,  a  mapping  will  be  provided  that  will  show  regions  of  agreement  and 
regions  of  conflict  for  two  different  methods  used  to  predict  handling  qualities.  Since  every 
point  in  a  region  represents  a  different  aircraft  transfer  function,  a  variable  stability  aircraft 
can  be  configured  to  one  of  these  transfer  functions.  Table  1.2  shows  the  operating  envelopes 
of  some  parameters  for  the  F-16  VISTA  (18).  A  point  in  the  region  of  conflict  could  be  ‘dialed 
in’  to  the  flight  control  system.  A  test  pilot  could  then  fly  the  simulated  aircraft  and  determine 
what  the  real  handling  qualities  level  is.  If  enough  cases  are  flight  tested  and  the  results  show 
a  trend,  a  decision  could  be  made  on  whether  one  of  the  methods  used  to  predict  handling 
qualities  needs  to  be  modified  or  eliminated.  This  thesis  lays  the  foundation  for  this  type  of 
future  analysis. 


Table  1.2  Operating  Envelopes  for  the  F-16  VISTA 


Parameter 

Range 

Short  Period  Natural  Frequency 

0  to  12  rad/sec 

Short  Period  Damping 

-0.1  to  1.1 

Aircraft  Time  Delay 

0.01  to  0.5  sec 

Experience  with  aircraft  operations  indicates  that  certain  flight  phases  require  more 
stringent  values  of  flying  qualities  parameters  (4:80).  A  given  flight  phase  can  have  an 
aircraft  normal  state  associated  with  it.  For  example,  the  flaps  and  gear  are  down  for  landing 
approaches  and  up  for  cruising  flight.  MIL-STD-1797A  defines  three  categories  of  flight 
phases  as  outlined  in  Table  1.3. 
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Table  1.3  Flight  Phase  Categories 


Category 

Description 

Example 

A 

nonterminal  flight  phases  that  require 
rapid  maneuvering,  precision  tracking, 
or  precise  flight-path  control 

air-to-air 

combat 

B 

nonterminal  flight  phases  that  require 
gradual  maneuvering  without 
precision  tracking 

climb 

C 

terminal  flight  phases  that  require 
gradual  maneuvering  with 
accurate  flight-path  control 

landing 

This  research  only  examines  the  Category  C  flight  phase.  However,  the  approach 
derived  in  Chapter  III  can  handle  any  flight  phase  category.  With  aircraft  classification  and 
flight  phase  known,  handling  qualities  levels  can  now  be  addressed. 

MIL-STD-1797A  defines  three  levels  of  handling  qualities  which  are  outlined  in  Ta¬ 
ble  1.4  (4:85).  These  definitions  are  based  on  the  Cooper-Harper  (2)  rating  scale  shown  in 
Fig.  1.4. 


Table  1.4  Handling  Qualities  Levels 


Level 

Meaning 

Description 

1 

satisfactory 

Flying  qualities  clearly  adequate  for  the  mission 
flight  phase.  Desired  performance  is  achievable 
with  no  more  than  minimal  pilot  compensation. 

2 

acceptable 

Flying  qualities  adequate  to  accomplish  the  mission 
flight  phase,  but  some  increase  in  pilot  workload 
and/or  degradation  in  mission  effectiveness  exists. 

3 

controllable 

Flying  qualities  are  such  that  the  aircraft  can  be 
controlled  in  the  context  of  the  mission  flight 
phase,  even  though  pilot  workload  is  excessive 
and/or  mission  effectiveness  is  inadequate. 

It  is  not  uncommon  to  have  feedback  control  systems  of  twentieth  order  or  more.  Writing 
a  specification  for  such  a  system  presents  a  problem.  Considerable  research  has  been  devoted 
to  reducing  the  order  of  these  high-order  feedback  control  systems  by  matching  frequency 
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Figure  1.4  Cooper-Harper  Scale 


responses  to  obtain  lower-order  equivalent  systems  (4:175).  Using  lower-order  equivalent 
systems  allows  the  application  of  well-established  boundaries  generated  by  classical  airplane 
data  to  be  extended  to  many  high  order  systems  (4:175).  With  respect  to  pitch  axis  control 
requirements,  boundaries  have  been  established  on  the  classical  modal  parameters  in  Eq.  (1.1). 

M  =  I<e  (Tei  ■  s  +  1)  (Tq2  ■  s  +  l)  e~^s 

fi(s)  (s2  +  2(pUp  'S-fWp)  (s2  +  2 CspLOsp  •  5  +  w2p) 


where, 

I<g  =  aircraft  gain 

Tgx  -  low-frequency  pitch  attitude  zero 

Tq2  =  high-frequency  pitch  attitude  zero 

tq  =  aircraft  time  delay 

CP  =  damping  ratio  of  the  phugoid  mode 

up  =  undamped  natural  frequency  of  the  phugoid  mode 

(sp  =  damping  ratio  of  the  short  period  mode 

usp  =  undamped  natural  frequency  of  the  short  period  mode 

s  =  Laplace  variable 


Equation  (1.1)  is  a  linearized,  reduced-order  model  of  the  actual  aircraft  response.  The  flight 
control  system  is  also  embedded  in  Eq.  (1.1).  In  most  cases,  the  phugoid  and  short-period 
modes  are  sufficiently  separated  that  further  order  reduction  is  possible.  Equation  (1.2), 
universally  recognized  as  the  pitch  model  of  short-period  dynamics,  may  normally  be  used  in 
place  of  Eq.  (1.1)  (4:175). 

m  =  Kt(Tt2  ■  s  +  (U) 

S(s)  a  (s2  +  2Cr“v  '  a  +  w lv ) 

For  example,  the  following  pitch  transfer  function  is  for  an  A-4  aircraft  with  a  feel 
system  and  actuators  flying  at  Mach  0.85  (6:9). 

d(s )  1.07745  •  104(s  +  2.085)(s  +  0.0287) 

8(s)  ~  (s2  +  36.4s  +  676)(s2  +  6.39s  +  54.02)(s2  +  0.00313s  +  0.00484) 
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Equation  (1.4)  is  the  short-period  approximation  of  Eq.  (1.3). 

0(s)  16.26(s  +  2.085)e-°-O9471  s 

6(4  _  s{s2  +  6.98s  +  51.29) 

Figure  1.5  compares  Bode  plots  of  the  original  higher  order  system  to  the  short-period  ap¬ 
proximation.  For  this  example,  the  short-period  approximation  describes  the  aircraft  pitch 
response  fairly  well  for  the  frequency  range  given.  Both  the  u )spTg2,(sp,  Tg  and  Neal-Smith 
criteria  are  only  valid  for  frequencies  up  to  10  radians  per  second  (4:190,237). 

12  Objectives 

The  overall  goal  of  this  thesis  is  to  compare  the  ojSpTg2  ,  (sp,  Tg  and  Neal-Smith  criteria. 
The  following  objectives  are  needed  to  accomplish  this  goal.  First,  a  computer  program 
that  calculates  the  pilot  model  parameters  for  the  Neal-Smith  criteria  needs  to  be  developed. 
Second,  a  computer  program  that  maps  areas  from  the  tospTg2 ,  (sp,  Tg  criteria  into  the  Neal- 
Smith  criteria  needs  to  be  developed.  Third,  regions  of  agreement  and  regions  of  conflict  need 
to  be  identified  on  the  Neal-Smith  criteria.  Forth,  a  computer  program  that  maps  these  regions 
back  into  the  u3pTg2,(sp,  t6  criteria  needs  to  be  developed.  Finally,  a  parametric  study  to 
identify  trends  needs  to  be  performed. 

1 .3  Assumptions 

The  first  assumption  is  the  aircraft  transfer  function  for  pitch  control  can  be  modeled 
by  a  short-period  approximation  which  includes  the  flight  control  system.  This  assumption 
is  only  needed  when  the  c ospTg2,(sp,  t8  criteria  is  mapped  into  the  Neal-Smith  criteria.  The 
algorithm  developed  in  Chapter  IB  for  determining  the  Neal-Smith  criteria  can  handle  any 
aircraft  transfer  function. 

The  second  assumption  is  the  pilot  can  be  modeled  by  a  variable  gain,  a  constant  time 
delay,  and  a  variable  first-order  compensation  network.  This  assumption  is  needed  for  the 
Neal-Smith  criteria. 
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Amplitude  (dB) 


1 .4  General  Approach 

The  lospTo2,(sp,  Te  criteria  involves  four  parameters.  Holding  Tq2  and  tq  fixed  while 
varying  cosp  and  (sp  forms  a  two-dimensional  region  that  can  be  mapped  into  the  Neal-Smith 
criteria.  For  example,  the  region  might  be  all  points  corresponding  to  level  1  handling  qualities. 
The  actual  values  of  all  four  parameters  completely  define  the  aircraft  transfer  function.  The 
remainder  of  the  problem  involves  solving  for  the  parameters  of  the  pilot  transfer  function. 

The  assumed  pilot  transfer  function  has  three  unknowns.  The  resulting  equations  are 
complicated  functions  of  frequency  with  no  closed-form  solution.  The  solution  is  found  by 
posing  the  problem  as  a  constrained  optimization  problem.  Once  all  three  unknowns  are 
determined,  the  mapping  of  the  original  point  used  to  define  the  aircraft  transfer  function  can 
be  accomplished.  The  procedure  is  repeated  for  any  point  contained  by  the  two-dimensional 
region  described  above.  Ultimately,  a  grid  of  points  representing  the  region  will  be  mapped. 

The  inverse  mapping  would  be  even  more  difficult.  A  point  from  the  Neal-Smith  criteria 
represents  the  pilot-aircraft  system  in  the  form  of  one  transfer  function.  The  individual  pilot 
and  aircraft  parameters  are  not  easily  distinguished.  Rather  than  trying  to  accomplish  this 
for  a  grid  of  points,  bookkeeping  of  the  original  mapping  is  used.  When  a  point  from  the 
ujspTg2 ,  (sp,  Tg  criteria  is  mapped  into  the  Neal-Smith  criteria,  information  on  the  Neal-Smith 
handling  qualities  level  is  stored  and  identified  with  that  point.  After  the  entire  grid  is  mapped, 
every  point  from  uspTg2,(sp,  Tg  is  assigned  a  handling  qualities  level  from  the  Neal-Smith 
criteria.  The  same  method  can  be  applied  to  the  Pilot-in-the-Loop  criteria. 

1.5  Sequence  of  Presentation 

Chapter  II  contains  an  in-depth  look  at  the  cospTg2 ,  (sp,  Tg,  Neal-Smith  and  Pilot-in-the- 
Loop  criteria.  Chapter  in  contains  a  detailed  description  of  the  approach  used  to  perform  the 
mappings.  Chapter  IV  contains  the  results  of  the  mappings  of  the  Category  C  flight  phase  of 
the  Class  II  and  Class  IV  aircraft.  Chapter  V  contains  conclusions  and  recommendations  for 
future  research. 
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II.  Pitch  Response  Criteria 


2.1  The  LospTo2,(sP,Te  Criteria 

Physically,  lospTo2  represents  the  lag  in  phase  (at  tosp)  or  time  between  aircraft  responses 
in  pitch  attitude  and  path.  If  Tg2  is  too  small  with  respect  to  usp,  the  path  and  attitude  response 
may  not  be  separated  enough  to  give  a  pilot  the  additional  cues  he  or  she  needs  to  control 
the  slower  path  loop  (4:192).  ujspTg2,  in  combination  with  (sp,  defines  the  shape  of  the 
attitude  frequency  response  (4:192).  If  the  aircraft  transfer  function  in  Eq.  (1.2)  is  known,  the 
handling  qualities  level  can  be  read  directly  from  Fig.  2.1  for  a  Category  C  flight  phase.  There 
is  no  upper  limit  placed  on  LospTg2.  However,  MIL-STD-1797A  only  plots  coSpTe2  UP t0  ten- 
Therefore,  the  highest  value  of  uspTg2  in  this  thesis  is  ten. 


Figure  2. 1  Handling  Qualities  Boundaries  for  the  u>spTe2 ,  (sp,  re  criteria 
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Table  2.1  Limits  on  rs 


LEVEL 

Tdmax  (seconds) 

1 

2 

3 

0.30 

Table  2.2  Minimum  Values  for  tosp  and  Maximum  values  for  Tq2 


LEVEL 

CLASS 

{yJsp)min 

(Te2  )max 

1 

I,  iv 

■Mi 

II,  III 

2 

i,  IV 

0.60 

■SB-'- 

II,  III 

7.14 

0.14 

2.1.1  Maximum  Values  for  re.  In  addition  to  satisfying  the  boundaries  in  Fig.  2.1, 
the  aircraft  equivalent  time  delay  (to)  must  not  be  more  than  a  certain  value.  MIL-STD-1797A 
maximum  values  for  t0  with  respect  to  handling  qualities  level  are  listed  in  Table  2.1. 


2.1.2  Minimum  Values  for  cosp  and  Maximum  values  for  Tg2.  MIL-STD-1797A 
places  limits  on  usp  and  Tg2  as  outlined  in  Table  2.2.  These  limits  are  dependent  on  aircraft 
classification  as  well  as  handling  qualities  level.  Normal  acceleration  response  to  attitude 
changes  is  a  primary  factor  affecting  the  pilot’s  perception  of  the  minimum  allowable  usp 
(4:192).  The  maximum  values  for  T0,2  are  based  on  an  approach  speed  of  135  kt  (4:206). 


2.2  The  Neal-Smith  Criteria 

The  Neal-Smith  criteria  is  not  as  straight  forward  as  the  uspTg2:(sp,  to  criteria.  The 
main  difference  is  the  Neal-Smith  criteria  uses  a  closed-loop  system  that  incorporates  a  pilot 
model.  Figure  2.2  shows  the  block  diagram  of  the  pilot-aircraft  system.  The  system  is  a 
single-input-single-output  system  that  uses  negative  unity  feedback.  Understanding  of  the 
following  terminology  is  needed  for  the  remaining  sections.  Yp  is  the  pilot  transfer  function. 
|  is  the  open-loop  transfer  function  of  the  airplane  plus  control  system,  Yc.  Equation  1.2  will 
be  used  for  this.  £  is  the  open-loop  transfer  function  of  the  airplane  plus  control  system  plus 
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Figure  2.2  Block  Diagram  of  the  Pilot-Aircraft  System  used  in  the  Neal-Smith  Criteria 
pilot  which  is  given  the  symbol  G. 

G  =  YPYC  (2.1) 


f-  is  the  closed-loop  transfer  function  of  the  airplane,  control  system,  and  pilot.  Equation  2.2 
represents  this  transfer  function  which  is  given  the  symbol  T. 


G 

1  +G 


(2.2) 


T  is  also  known  as  the  complimentary  sensitivity  (7).  The  magnitude  of  any  transfer  function 
will  be  designated  |  ■  |,  and  the  phase  of  any  transfer  function  will  be  designated  L  (•). 

2.2.1  A  Pilot’s  View  of  Good  Tracking  Performance.  What  is  a  pilot  trying  to 
accomplish  when  he  or  she  adapts  to  an  aircraft  configuration?  Neal  and  Smith  examined  pilot 
comments  and  led  to  the  following  conclusion.  The  pilot  wants  to  acquire  the  target  quickly 
and  predictably,  with  minimum  of  overshoot  and  oscillation  (12:39).  References  (9)  and  (11) 
place  the  following  mathematical  relationships  on  pilot  comments: 

( 1 )  The  pilot  tries  to  achieve  a  particular  value  of  the  open-loop  gain  crossover  frequency, 
ujc  (the  frequency  at  which  £  is  0  dB). 

(2)  The  pilot  tries  to  minimize  any  closed-loop  ‘droop’  (hold  £•  near  0  dB,  for  u;  <  loc). 

(3)  The  pilot  tries  to  maintain  good  high-frequency  stability  by  keeping  the  damping 
ratio  of  any  closed-loop  oscillatory  modes  greater  than  0.35,  and  by  maintaining  a  phase 
margin  of  60  to  1 10  degrees. 


20 


Figure  2.3  Designation  of  Droop  and  Resonance 

Neal  and  Smith  substituted  another  closed-loop  relation  in  place  of  (1).  First  a  closed- 
loop  parameter  similar  to  the  crossover  frequency,  u>c,  needed  to  be  defined.  Neal  and  Smith 
defined  a  parameter  called  closed-loop  bandwidth  (oobw)  to  be  the  frequency  at  which  the 
closed-loop  phase  (/£)  is  -90  degrees  (12:41). 

Q 

L—  =  —90°  for  to  =  u>bw  (2.3) 


Neal  and  Smith  determined  that  the  pilot  wants  to  minimize  closed-loop  droop  for  frequencies 
below  u>bw-  For  frequencies  above  cobw,  the  pilot  wants  to  minimize  the  resonant  peak 
(12:42).  Figure  2.3  illustrates  droop  and  resonance. 

2.2.2  Tracking  Performance  Standards.  Neal  and  Smith  defined  the  minimum 
value  for  closed-loop  droop  to  be  -3  dB  (12:43). 


>  —3 dB  for  to  <  cobw 


(2.4) 
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Figure  2.4  illustrates  this  constraint  along  with  closed-loop  bandwidth. 

Another  way  to  look  at  these  constraints  with  regards  to  open-loop  information  is  via  a 
Nichol’s  chart.  A  Nichol’s  chart  shows  closed-loop  boundaries  overlaid  on  a  plot  of  open-loop 
magnitude  versus  open-loop  phase.  Figure  2.5  illustrates  the  closed-loop  bandwidth  and  droop 
requirements  on  aNichol’s  Chart.  Also  shown  on  Fig.  (2.5)  are  the  upper  boundaries  according 
to  MIL-STD-1797A.  These  upper  boundaries  will  be  discussed  further  in  the  section  on  the 
Pilot-in-the-Loop  criteria.  So  far,  only  the  aircraft  transfer  function  has  been  discussed.  The 
pilot  transfer  function  will  be  the  focus  of  the  next  section. 

2.2.3  The  Pilot  Model.  Neal  and  Smith  require  the  pilot  transfer  function,  Yp,  to  be 

=  A'P(^^  +  l)p'  (2.5) 

The  pilot  includes  a  variable  gain  (Kp),  a  constant  time  delay  (rp),  and  a  variable  first-order 
compensation  network  involving  the  parameters  TPl  and  TP2.  Past  work  suggests  tp  may  vary 
with  airplane  dynamics  and  will  usually  lie  between  0.2  and  0.4  seconds  (12:38).  Neal  and 
Smith  assumed  tp  to  be  fixed  at  0.3  seconds  (12:38). 

tp  =  0.3  ( seconds )  (2.6) 

This  assumption  will  be  used  for  this  thesis.  The  remaining  parameters  I<p,  TPl ,  and  TP2  will 
need  to  be  solved.  So  far,  no  guidelines  or  criteria  have  been  given  for  resonance  limits.  The 
closed-loop  droop  must  be  greater  than  -3  dB  and  the  closed-loop  phase  must  be  —90°  when 
to  =  ujbw-  The  next  section  will  examine  closed-loop  resonance  and  pilot  phase  contribution. 

2.2.4  Neal- Smith  Boundaries  for  Handling  Qualities.  Neal  and  Smith  examined 
a  great  deal  of  pilot  comments  and  grouped  the  comments  as  a  function  of  closed-loop  and 
open-loop  parameters.  Figure  2.6  is  the  result  of  this  examination  (12:94).  Closed-loop 
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Figure  2.5  Nichol’s  Chart  with  MIL-STD-1797A  Boundaries 
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Figure  2.6  Summary  of  Pilot  Comments 


resonance  is  plotted  against  the  pilot  phase  angle,  which  is  defined  as 


Pilot  Phase  Angle  =  L 


( '  Tp\  +  1 
\ju>  ■  TP2  +  1 


(2.7) 


From  Fig.  2.6,  Neal  and  Smith  proposed  the  boundaries  shown  in  Fig.  2.7.  These  boundaries, 
along  with  closed-loop  droop  and  bandwidth,  will  serve  as  the  design  criteria  when  solving  for 
Kp,  TPl ,  and  TP2.  A  unique  solution  is  desired  for  mapping.  The  remainder  of  Chapter  II  and 
the  beginning  of  Chapter  III  will  show  how  to  reduce  this  three-degree-of-freedom  problem 
into  a  one-dimensional  search. 


2.2.5  Lead  Compensation.  The  reader  is  referred  to  (7)  for  additional  information 
on  various  forms  of  compensation.  The  following  transfer  function  can  be  used  to  explain 
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Figure  2.7  Neal-Smith  Boundaries 


both  lead  and  lag  compensation. 


D(s) 


Kg{Td  -3  +  1) 

(aTd  -3  +  1) 


(2.8) 


For  lead  compensation,  a  <  1.  Equation  (2.8)  is  similar  to  the  variable  first  order  com¬ 
pensation  network  used  for  the  pilot  model  in  Eq.  (2.5).  Lead  compensation  adds  phase  to 
the  system  at  frequencies  above  the  breakpoint  1/7).  At  the  same  time,  the  magnitude  of 
the  compensation  continues  to  grow  with  frequency.  For  nonzero  a,  the  magnitude  of  the 
compensation  levels  off  at  frequencies  above  l /aTd.  In  relation  to  the  pilot  model,  a  <  1 
implies  TPl  >  TP2. 

The  effects  of  lead  compensation  can  also  be  seen  on  a  Nichol’s  chart.  Curve  A  in 
Fig.  2.8  meets  the  closed  loop  bandwidth  requirement,  but  the  resonance  is  greater  than  9  dB. 
Lead  compensation  causes  the  lower  part  of  curve  A  to  shift  upward  and  to  the  right  (shaped 
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more  like  curve  C)  (12:47).  The  pilot  gain  can  then  be  adjusted  to  shift  curve  A  downward  to 
again  meet  the  closed-loop  bandwidth  requirement,  while  at  the  same  time  reducing  resonance. 

A  ‘target  zone’  can  be  formed  from  the  closed-loop  bandwidth  and  droop  require¬ 
ments.  The  —90°  closed-loop  phase  line  on  Fig.  2.8  terminates  at  the  point  where  open-loop 
phase=— 180°.  The  —90°  closed-loop  phase  line  and  the  -3  dB  droop  line  intersect  at  a  point 
where  the  open-loop  phase=— 125.3°.  As  long  as  the  point  where  a;  =  ujbw  satisfies  —180°  < 
open-loop  phase  <  -125.3°,  the  closed-loop  bandwidth  requirement  can  be  met  by  adjusting 
the  gain.  If  the  point  where  u  =  toBw  corresponds  to  a  point  where  the  non-compensated 
open-loop  phase  <  —180°,  lead  compensation  is  needed  to  shift  that  point  to  the  right.  The 
point  would  also  shift  up.  The  gain  would  need  to  be  adjusted  to  bring  the  point  back  down 
to  the  original  magnitude.  The  idea  is  to  tune  the  parameters  to  give  more  movement  to  the 
right  than  up.  The  net  result  would  be  an  overall  shift  to  the  right.  Another  time  when  lead 
compensation  is  needed  occurs  when  the  closed-loop  bandwidth  requirement  is  met  but  the 
droop  is  above  the  -3  dB  droop  line.  This  was  already  shown  for  curve  A  in  Fig.  2.8.  If  the 
point  where  to  =  lobw  is  to  the  right  of  the  target  zone,  lag  compensation  will  be  needed. 

Before  discussing  lag  compensation,  the  calculation  of  non-compensated  open-loop 
phase  needs  to  be  addressed.  IYC  is  the  open-loop  phase  of  the  aircraft  with  the  aircraft  time 
delay  (t#)  included.  The  pilot  also  has  a  time  delay.  Since  the  two  time  delays  are  multiplied 
when  forming  G,  they  can  be  replaced  by  a  combined  time  delay,  r,  as  in  Eq.  (2.9). 

e-res  .  e-rp-s  _  e-(re+r _  e-r«  (2.9) 

where  r  =  rg  +  tp 

When  calculating  the  non-compensated  open-loop  phase,  the  combined  time  delay  will  be 
used.  Rather  than  introduce  a  new  symbol,  LYC  will  always  include  the  phase  contributed  by 
e~Tp's. 
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o  DESIGNATES  (0  =  eoBW 


2.2.6  Lag  Compensation.  For  lag  compensation,  a  >  1  or  TPl  <  TV2.  Lag 
compensation  does  the  opposite  of  what  lead  compensation  does.  Phase  is  subtracted  and 
magnitude  decreases.  This  can  also  be  seen  on  aNichol’s  chart.  Curve  B  in  Fig.  2.8  meets  the 
closed-loop  droop  requirement,  but  it  does  not  meet  the  closed-loop  bandwidth  requirement. 
Lag  compensation  causes  the  lower  part  of  curve  B  to  shift  downward  and  to  the  left  (shaped 
more  like  curve  C)  (12:49).  The  pilot  gain  can  be  adjusted  to  shift  curve  B  downward  to 
meet  the  closed-loop  bandwidth  requirement,  reduce  resonance,  and  still  meet  the  closed-loop 
droop  requirement.  In  addition  to  being  to  the  right  of  the  target  zone  when  <x>  =  lobw,  lag 
compensation  is  also  needed  when  closed-loop  bandwidth  is  met  but  the  droop  is  below  the 
-3  dB  droop  line.  Deciding  whether  to  use  lead  or  lag  compensation  is  an  important  part  of 
the  approach  in  Chapter  III. 

Neal  and  Smith  state  that  the  lowest  resonance  occurs  when  the  closed-loop  droop 
equals  -3  dB  while  meeting  the  closed-loop  bandwidth  requirement  simultaneously  (12:47- 
49).  Curve  C  in  Fig.  2.8  illustrates  their  statement.  Adding  lead  compensation  to  curve  C 
would  cause  the  resonance  to  increase  if  the  closed-loop  droop  is  held  at  -3  dB.  Adding  lag 
compensation  to  curve  C  would  cause  the  resonance  to  increase  if  the  closed-loop  bandwidth 
is  held  at  u>Bw-  Changing  Eq.  (2.4)  to  an  equality  constraint  provides  a  unique  solution 
for  the  pilot  model.  A  unique  solution  is  needed  to  perform  a  mapping.  Searching  over  a 
three-dimensional  design  space  for  a  unique  solution  is  not  practical.  Reducing  the  order  of 
the  problem  is  the  next  step. 

2.2.7  ‘Optimum’  Pilot  Compensation.  Neal  and  Smith  define  ‘optimum’  lead 
compensation  as  providing  the  most  positive  increase  in  Lj-  for  the  least  flattening  of  the 
amplitude-phase  curve,  in  the  general  vicinity  of  ujBw  (12:50).  The  following  analysis  further 
explains  the  meaning  of  the  last  statement.  Equation  (2.10)  represents  the  amplitude  (in  dB) 
of  the  pilot  model  with  Kp  —  1.  Equation  (2.11)  represents  the  phase  (in  radians)  of  the  pilot 
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Figure  2.9  Amplitude-Phase  Curve  for  Lead  Compensation 


model. 


20  log 


ju  ■  TVl  +  1 
ju  •  TP2  +  1 


10  •  log  (u2T2pi  +  l)  -  10  •  log  (lo2T22  +  l) 


(dB)  (2.10) 


fju-TPl  +  1\ 
\ju  ■  TP2  +  1 ) 


tan  1  (uTp  j )  —  tan 


-l 


LJPi 


(2.11) 


The  flattening  tendency  is  related  to  the  fact  that  the  increment  in  open-loop  amplitude 
contributed  by  lead  compensation  is  positive  and  increases  with  to.  Figure  2.10  illustrates  this 
for  various  values  of  TP2  /  TPl .  Least  flattening  can  be  thought  of  as  providing  the  most  amount 
of  phase  for  a  given  amplitude.  Figure  2.9  shows  that  for  any  amplitude,  the  maximum  amount 
of  phase  occurs  when  TP2/TPl  =  0.  Another  way  to  explain  flattening  is  by  examining  the 
slopes  in  Fig.  2.10.  Equation  2.12  is  the  derivative  of  Eq.  (2.10)  with  respect  to  logo;. 
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Figure  2.10  Magnitude  Contribution  of  Lead  Compensation 
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Figure  2.11  Phase  Contribution  of  Lead  Compensation 
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Figure  2.12  Magnitude  Slope  and  Phase  Contributed  by  Lead  Compensation 


Slope  is  given  the  symbol  m. 


m  = 


d  (log  w) 
20 


20  log 


ju>  •  Tp j  +  1 


juj  ■  TP2  +  1 

^(l+w-’r-2) 


1  - 


(l+o;-2Tpf)  [  (rp22r-2  +  w~2T~2) 


(dB/decade)  (2.12) 


Figure  2.12  plots  this  slope  against  phase.  Least  flattening  for  a  given  amount  of  phase 
corresponds  to  the  least  positive  value  in  slope.  Referring  again  to  Fig.  2.10,  the  slopes 
increase  to  maximum  values  at  some  intermediate  frequency,  then  decrease.  To  determine  this 
extremum,  the  partial  derivative  of  Eq.  (2.11)  with  respect  to  co  is  set  equal  to  zero. 


A  [/  (iu'TPi  +  1\ 
du  [  \ju  -TP2  +  l) 


T  - 

Jp\ 

f rp  rp2  rp2  rp  \ 

k1V11V2  -LP1-LP2J 

1 

(1+^) 

1  1 

(l+^J 

1 

(2.13) 
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This  corresponds  to  setting  the  numerator  in  Eq.  (2.13)  to  zero. 


TP  j  —  TP2  +  to2 


=  0 


(2.14) 


The  extremum,  u)ext,  becomes 


M ext  —  i 


\JTP1  '  TV2 


(2.15) 


To  determine  whether  uext  is  a  maximum  or  minimum,  the  second  partial  derivative  of 
Eq.  (2.11)  is  taken. 


duj 2 


'ju  '  TPl  +  1 
Jw  •  Tpi  + 1 , 


2u;Tp2 


2^ 


(l  +“2TP\)  (l  +^)' 


(2.16) 


Setting  u>  =  uext  gives 


[  \iw  ■  TP2  +  1/ 


2 1%  (l  +  TnT-')2  -  2TI  (l  +  T-'TKf 

(l  +  TnT~')2  (l  +  T-'T„Y 


(2.17) 


Since  the  denominator  in  Eq.  (2.17)  is  positive,  the  numerator  determines  the  sign.  If  Eq.  (2.17) 
is  positive,  u>ext  is  a  minimum.  If  Eq.  (2.17)  is  negative,  u)ext  is  a  maximum.  For  lead 
compensation,  0  <  TP2/TPl  <  1.  Figure  2.13  illustrates  the  sign  of  the  numerator  as  a 
function  of  TV2/TVl.  The  graph  shows  that  Eq.  (2.17)  is  negative  in  the  region  for  lead 
compensation.  Thus,  uext  is  the  maximum  for  lead  compensation.  This  can  be  seen  on 
Fig.  2.12.  Slope  and  phase  increase  until  reaching  the  center  frequency,  then  decrease  along 
the  same  path  back  to  zero.  While  coext  forms  an  upper  bound  to  Fig.  2.12,  it  is  not  the  solution 
for  ‘optimum’  lead  compensation  as  defined  by  Neal  and  Smith.  Instead,  the  most  positive 
phase  for  a  given  positive  slope  is  always  when  TP2/TPl  =  0.  Thus,  the  ‘optimum’  lead 
compensation  is  pure  lead,  or 

TP2  =  0.  (2.18) 
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Figure  2.13  2nd  partial  of  phase  as  a  function  of  TP2  /  Tn 

Neal  and  Smith  define  ‘optimum’  lag  compensation  as  providing  the  most  steepening  of 
the  amplitude  phase  curve  for  the  least  negative  increase  in  L  in  the  vicinity  of  uBw  (12:50). 
For  lag  compensation,  TP2/TPl  >  1.  An  analysis  similar  to  the  one  for  lead  compensation 
helps  to  explain  the  meaning  of  Neal  and  Smith’s  definition.  The  steepening  tendency  is 
related  to  the  fact  that  the  increment  in  open-loop  magnitude  contributed  by  lag  compensation 
is  negative  and  decreases  with  u>.  Figure  2.15  illustrates  this  for  various  values  of  TP2/TPl. 
The  steepest  point  on  each  curve  occurs  at  u>ext .  In  terms  of  slope,  ‘optimum’  lag  compensation 
provides  the  most  negative  value  of  m  for  the  least  negative  value  of  phase  as  seen  in  Fig.  2.17. 
Slope  and  phase  decrease  to  minimum  values  at  toext,  then  increase  along  the  same  path  back 
to  zero.  Referring  gain  to  Fig.  2.13,  Eq.  (2.17)  is  positive  when  TP2/TPl  >  1.  Therefore, 
uext  is  the  minimum  for  lag  compensation.  Unlike  lead  compensation,  u>ext  is  the  solution  for 
‘optimum’  lag  compensation.  The  most  negative  slope  for  a  given  negative  phase  is  always 
at  the  center  frequency.  Since  the  primary  area  of  interest  is  in  the  vicinity  of  c oBw,  Neal  and 
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Figure  2.14  Amplitude-Phase  Curve  for  Lag  Compensation 


Smith  set  u>ext  =  ojbw-  Thus,  ‘optimum’  lag  compensation  occurs  when 


U>BW  =  ± 


\jTVl  •  Tp 2 


(2.19) 


Solving  Eq.  (2.19)  for  TP2  yields 


T  —  — 

P2  r, 


Pi  '  UBW 


(2.20) 


If  either  ‘  optimum  ’  lead  compensation  or  ‘optimum’  lag  compensation  is  used,  Eq.  (2. 1 8) 
or  Eq.  (2.20)  reduce  the  order  of  the  three-dimensional  problem  down  to  two. 
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Figure  2.17  Magnitude  Slope  and  Phase  Contributed  by  Lag  Compensation 


2.3  The  Pilot-in-the-Loop  Criteria 


The  Pilot-in-the-Loop  criteria  is  a  modified  version  of  the  Neal-Smith  criteria.  The 
differences  are  outlined  in  this  section.  While  Neal  and  Smith  require  only  one  form  of  pilot 
model,  MIL-STD-1797A  offers  two  pilot  models  for  use  with  the  Pilot-in-the-Loop  criteria. 
The  first  pilot  model  is  the  same  one  Neal  and  Smith  use  in  Eq.  (2.5).  The  second  pilot  model 
is  used  when  there  is  no  free  s  in  the  denominator  of  the  aircraft  transfer  function  (4:237).  A 
free  s  is  desired  in  G  to  make  the  system  a  type  I  system.  Type  I  systems  have  zero  steady-state 
error  to  step  inputs  (7:108).  Equation  (2.21)  can  be  used  with  the  Pilot-in-the-Loop  criteria 
when  Yc  has  no  free  s. 


Kp(5s  +  l)(rpi  •  3  +  l)e~^s 
s  '  (Tp2  •  s  +  1) 


(2.21) 


Since  the  short  period  approximation  will  be  used  for  Yc,  Eq.  (2.5)  will  be  used  for  Yp. 

The  second  difference  is  in  the  value  of  the  pilot  time  delay,  tp.  Neal  and  Smith  used 
tp  =  0.3  seconds.  MIL-STD-1797A  fixes  tp  at  0.25  seconds  (4:237).  This  difference  can  be 
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ignored  if  the  pilot  and  aircraft  time  delays  are  combined.  Data  generated  for  one  value  of  r  can 
be  used  to  analyze  both  the  Neal-Smith  criteria  and  the  Pilot-in-the-Loop  criteria.  The  aircraft 
time  delay,  however,  would  be  different  for  each  criteria.  For  example  if  r  =  0.5  seconds, 
to  would  be  0.25  seconds  for  Pilot-in-the-Loop  analysis  and  0.20  seconds  for  Neal-Smith 
analysis. 

The  final  difference  between  the  Neal-Smith  and  Pilot-in-the-Loop  criteria  is  limits 
on  closed-loop  resonance.  MIL-STD-1797A  puts  the  following  limits  on  the  resonant  peak 
values. 

A  bandwidth,  defined  by  a  closed-loop  phase  of  -90  degrees,  of  2.5  rad/sec 
[Category  C]  shall  be  attainable  with  closed-loop  droop  no  more  than  -3  dB  for 
Levels  1  and  2  and  closed-loop  resonance  no  greater  than  3  dB  for  Level  1, 9  dB 
for  Level  2  over  the  frequency  range  from  0  to  10  rad/sec. 

(4:237) 

These  limits  are  displayed  in  Table  2.3  and  are  less  restrictive  than  the  Neal-Smith  criteria.  In 
addition,  there  is  no  consideration  given  to  pilot  phase  angle. 


Table  2.3  Limits  on  Closed-Loop  Droop  and  Resonance 


LEVEL 

MAX  DROOP 
(dB) 

MAX  RES 
(dB) 

C-L  PHASE 

@U>BW 

1 

-3 

3 

-90° 

2 

-3 

9 

-90° 

2.4  Summary 

The  ojspTe2 ,  (sp,  Tg,  Neal-Smith,  and  Pilot-in-the-Loop  criteria  were  explained  in  detail. 
The  LospTg2 ,  (sp,  tq  criteria  depends  only  on  the  aircraft.  Limits  on  u>sp,  Tg2,  and  re  were  given. 
If  tq  and  To2  are  fixed,  the  parameters  u>sp  and  (sp  can  be  varied  to  form  a  two  dimensional 
region  for  mapping.  Every  point  in  the  region  defines  an  aircraft  transfer  function,  Yc. 
Once  Yc  is  known,  a  search  for  a  pilot  transfer  function  can  begin.  The  form  of  the  pilot 
model  has  three  unknown  parameters  (Kp,  TPl,  and  TPl).  The  Neal-Smith  criteria  and  the 
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Pilot-in-the-Loop  criteria  are  identical  with  respect  to  the  unknown  pilot  model  parameters. 
They  differ  in  the  value  of  pilot  time  delay  (rp).  The  Neal-Smith  criteria  fixes  tp  at  0.3 
seconds,  while  the  Pilot-in-the-Loop  criteria  fixes  tp  at  0.25  seconds.  Neal  and  Smith  also  use 
different  boundaries  to  distinguish  between  handling  qualities  levels.  Pilot  phase  angle  is  taken 
into  account,  whereas  the  Pilot-in-the-Loop  criteria  only  evaluates  the  maximum  closed-loop 
resonance.  Despite  these  differences,  the  same  approach  will  be  used  to  determine  the  pilot 
model  parameters.  Depending  on  whether  lead  or  lag  compensation  is  needed  determines  if 
Eq.  (2.18)  or  Eq.  (2.20)  is  used.  Either  equation  reduces  the  problem  to  a  two-dimensional 
problem  in  terms  of  Kp  and  TPl .  In  Chapter  III,  an  equation  will  be  derived  for  the  closed-loop 
bandwidth  requirement.  This  equation  will  serve  to  reduce  the  remaining  two-dimensional 
problem  down  to  a  one-dimensional  search. 


III.  Approach 


This  chapter  outlines  the  approach  used  to  map  the  ojSpTd2,(Sp,Te  criteria  into  the 
Neal-Smith  and  Pilot-in-the-Loop  criteria.  Neal  and  Smith  developed  a  graphical  method  for 
determining  the  pilot  model  parameters.  This  method  is  discussed  followed  by  an  alternate 
method  that  makes  use  of  numerical  optimization.  Finally,  a  step-by-step  procedure  is  given. 

3.1  Neal-Smith  Graphical  Method 

Computer  power  and  numerical  optimization  techniques  have  increased  since  Neal 
and  Smith  developed  their  criteria  in  1970  (12).  Nonetheless,  they  developed  a  method  for 
determining  the  pilot  model  parameters.  The  method  involved  iteration  based  on  graphical 
interpretations.  The  first  step  was  plotting  the  Nichol’s  chart  of  Yc  as  in  Fig.  2.8.  By  visual 
inspection,  they  decided  whether  lead  or  lag  compensation  was  needed.  An  initial  guess  at  the 
unknown  parameters  for  either  ‘optimal’  lead  compensation  or  ‘optimal’  lag  compensation 
was  made.  Another  Nichol’s  chart  was  generated,  this  time  of  the  pilot-aircraft  system  (G)  as 
in  Fig.  2.5.  The  pilot  gain,  I<p,  was  adjusted  to  meet  the  closed-loop  bandwidth  requirement. 
If  the  closed-loop  droop  was  greater  than  -3  dB,  additional  lead  was  needed  to  meet  the  closed- 
loop  droop  exactly.  If  the  closed-loop  droop  was  less  than  -3  dB,  additional  lag  was  needed. 
They  would  then  iterate  on  this  process  until  both  the  closed-loop  bandwidth  and  the  closed- 
loop  droop  requirements  were  met  simultaneously.  Each  iteration  required  looking  at  plots  in 
order  to  make  decisions  about  the  next  guess.  Although  they  derived  a  way  to  make  the  initial 
guess  close  to  the  solution,  the  process  still  involved  iterations  based  on  visual  examination  of 
plots.  Performing  this  process  for  a  grid  of  1000  different  aircraft  transfer  functions  would  be 
a  formidable  task.  Developing  a  method  that  does  not  require  human  interaction  at  each  step 
became  a  primary  goal  of  this  research.  The  remainder  of  this  chapter  outlines  an  alternate 
method  where  the  problem  is  posed  as  a  constrained  optimization  problem. 
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Table  3.1  General  Problem  Statement 


Minimize: 

/(*) 

x  6  3? 

objective  function 

Subject  to: 

9i(x)  <  0 

2  =  1,/ 

inequality  constraints 

hj(x)  =  0 

j  =  l,m 

equality  constraints 

x[<xk<  xuk 

k  =  1,  n 

side  constraints 

3.2  Sequential  Quadratic  Programming 

Since  there  are  criteria  (i.e.  constraints)  on  this  problem,  a  constrained  optimization 
routine  can  be  used  to  solve  for  the  unknown  pilot  model  parameters.  The  reader  is  referred 
to  Vanderplaats  (16)  for  various  methods  of  numerical  optimization.  The  author  chose  to 
use  Sequential  Quadratic  Programming  (SQP).  “SQP  methods  represent  state-of-the-art  in 
nonlinear  programming  methods”  (8:2-23).  Schittowski  (15)  has  implemented  and  tested  a 
version  of  SQP  that  out  performs  every  other  method  tested  in  terms  of  efficiency,  accuracy, 
and  percentage  of  successful  solutions  (over  a  large  number  of  test  problems).  The  remainder 
of  this  section  provides  a  brief  overview  of  SQP. 

The  general  problem  statement  is  shown  in  Table  3.1  where  *  is  the  design  vector,  / 
is  the  objective  function,  h  is  a  vector  of  equality  constraints,  and  g  is  a  vector  of  inequality 
constraints.  Elements  in  the  design  vector,  xk,  also  have  side  constraints.  The  principal  idea 
behind  SQP  is  the  formulation  of  the  quadratic  programming  (QP)  sub-problem  based  on  a 
quadratic  approximation  of  the  Lagrangian  function,  L,  given  by  Eq.  (3.1). 

I  m 

L{x,  A)  =  f(x)  +  xiSiix)  +  x(i+j)hAx)  (3-1) 

!=1  j  —  1 

The  symbol  A  represents  the  Lagrange  multiplier  of  each  constraint.  If  a  constraint  is  not  active 
(e.g.  g\  <  0),  the  corresponding  Lagrange  multiplier  must  be  zero.  The  QP  sub-problem 
is  obtained  by  linearizing  the  nonlinear  constraints.  The  QP  problem  statement  is  shown  in 
Table  3.2.  The  subscript  b  represents  iteration  number.  This  sub-problem  can  be  solved  using 
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Table  3 .2  QP  Problem  Statement 


Minimize:  Hbdb  +  S7f{xb)Tdb  x  e  3? 

Subject  to:  \jgi(x)Tdh  +  gi(x)  <0  i  =  1,  / 
S/hj(x)Tdb  +  /ij(aj)  =  0  j  =  l,m 


any  QP  algorithm.  The  solution,  d&,  is  used  to  form  a  new  iterate  *&+i  as  in  Eq.  (3.2). 


Xb- 1-1  —  Xb  T"  C*bdb 


(3.2) 


The  step  length,  ab,  is  determined  by  an  appropriate  line  search  procedure  such  that  a  sufficient 
decrease  in  a  merit  function  is  obtained.  The  merit  function  \P(a;)  in  Eq.(  3.3)  will  be  used. 

I  m 

$(«)  =  /(*)  +  Sr<max(°’&'(*))  +  J2r(i+j)hAx)  (3-3) 

*= 1  j= 1 


The  variable,  r,  is  a  penalty  parameter.  The  matrix  Hb  from  the  QP  sub-problem  is  a  positive 
definite  approximation  of  the  Hessian  matrix  of  the  Lagrangian  function  in  Eq.  (3.1).  Hb 
will  be  updated  using  the  Broyden,  Fletcher,  Goldfarb,  and  Shanno  (BFGS)  method  shown  in 
Eq.  (3.4). 


Hb+ 1  —  H  + 


ml 

<ibsb 


sTbHtsb 


where  sb  =  *&+i  —  xb 


96  =  V/(*6+i)  -  V/(*&) 


(3.4) 


3.3  ‘Optimizing’  the  Pilot  Model 

SQP  will  be  used  to  solve  for  the  pilot  model  parameters.  Strictly  speaking,  SQP  is 
not  ‘optimizing’  the  pilot  model.  Instead,  SQP  is  used  to  solve  the  complicated  equations 
used  to  determine  TPx  and  TP2  according  to  the  Neal-Smith  definitions  of  ‘optimum’.  Neal 
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and  Smith  state  that  the  lowest  resonance  occurs  when  the  closed-loop  droop  equals  -3  dB 
while  meeting  the  closed-loop  bandwidth  requirement  simultaneously  (12:47-49).  In  other 
words,  they  assume  the  lowest  resonance  always  occurs  when  the  droop  inequality  constraint 
is  active. 

Edkins  developed  a  computer  program  that  automated  the  original  Neal-Smith  method 
(6:14).  Given  a  range  on  TP1  and  TP2,  the  program  calculates  handling  qualities  levels  for 
various  combinations  in  the  search  interval.  Based  on  the  these  results,  the  ranges  on  TP1 
and  TP2  are  modified  until  a  satisfactory  result  is  obtained.  The  program  still  involves  human 
interaction  at  each  iteration  step.  Edkins  even  stated  that  “the  iteration  process  is  fairly 
lengthy”  (6:14). 

Since  the  pilot  wants  to  minimize  closed-loop  resonance  while  keeping  droop  above 
the  -3  dB  line,  the  problem  can  be  posed  in  the  general  format  of  Table  3.1.  In  other  words, 
minimize  closed-loop  resonance  subject  to  droop>-3  dB.  This  is  the  same  as  minimizing  the 
oo-norm  of  the  complimentary  sensitivity,  T.  If  the  goal  is  to  ‘optimize’  the  pilot  model,  why 
restrict  the  form  of  compensation  to  either  ‘optimal’  lead  or  ‘optimal’  lag  defined  by  Eq.  (2.18) 
and  Eq.  (2.20)?  There  could  be  a  combination  of  TPl  and  TP2  that  leads  to  an  overall  better 
solution.  The  problem  then  becomes  a  two-dimensional  search  for  the  best  combination  of 
Tp j  and  TP2  that  minimizes  closed-loop  resonance  subject  to  the  droop  inequality  constraint. 
This  approach  has  the  disadvantage  of  requiring  two  initial  guesses.  Since  SQP  looks  for 
the  minimum  solution  in  vicinity  of  the  initial  guesses,  there  is  also  a  greater  chance  of  SQP 
returning  a  local  minimum.  SQP  had  to  be  frequently  restarted  with  different  initial  guesses 
to  obtain  smaller  values  of  the  objective  function.  Edkins’  program  was  run  a  few  times  to 
determine  acceptable  initial  guesses.  This,  however,  greatly  increased  computational  time. 
Hence,  an  alternate  formulation  was  sought. 

Using  the  Neal-Smith  assumption  that  the  lowest  resonance  occurs  when  both  constraints 
are  met  with  equality  provides  a  unique  solution.  The  issue  of  whether  a  solution  is  a  local 
or  global  minimum  goes  away.  The  objective  is  to  drive  the  droop  to  -3  dB  exactly  while 
meeting  the  closed-loop  bandwidth  requirement.  The  remainder  of  this  chapter  outlines  this 
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type  of  approach.  Kp  and  TPl  are  still  unknown.  The  next  section  derives  an  equation  for  I<p 
which  depends  only  on  TP1. 

3.4  Closed-Loop  Bandwidth 

Figure  2.2  can  be  manipulated  to  look  like  Fig.  3.1  if  the  pilot  gain  ( Kp )  and  the  aircraft 
gain  (K$)  are  separated  from  both  transfer  functions  and  given  their  own  block  K. 


Figure  3.1  Block  Diagram  of  Pilot- Aircraft  System  with  Separate  Combined  Gain 


This  does  not  change  the  problem.  Equation  (3.5)  defines  the  new  combined  gain.  Equa¬ 
tion  (3.6)  becomes  the  pilot  transfer  function,  and  Eq.  (3.7)  becomes  the  aircraft  transfer 
function. 


K  =  Kp  ■  I<g 

Y  _  (^pi  ’  s  1) 

(TP2-s  +  1) 

(Tg2-s  +  l)e~^ 
s  (s2  -F  ‘2(>spu)sp  •  s  -f- 


(3.5) 

(3.6) 

(3.7) 


The  pilot-aircraft  system  still  is  given  the  symbol  G,  but  Eq.  (3.8)  is  used. 


G  =  K  •  YPYC 


(3.8) 


In  the  frequency  domain,  G  can  be  represented  by  its  magnitude  multiplied  by  its  phase. 

G  =  \G\  ■  cos  (j>G  +  j  \G\  ■  sin  <j>G  (3.9) 
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where  \G\=  Magnitude  of  G  and  <j)G  =  Phase  of  G 


The  closed-loop  transfer  function  becomes 


G  \G\  ■  cos  <f>G +j\G\- sin  <f>G 


1  +  G  l  +  \G\  ■  cos  (f>G  +  j  \G\  ■  sin  (j)G 


(3.10) 


Multiplying  the  numerator  and  denominator  of  Eq.  (3.10)  by  the  complex  conjugate  of  its 
denominator  yields 

T  _\G\-  cos  (f>G  +  \G\ 2  +  j  |G]  •  sin  <pG 
1+2  |G|  •  cos  (f)G  +  |G| 

The  phase,  or  argument,  of  T  becomes 


<j>CL  =  arg  [T]  =  arctan 


Imag(T ) 


Real(T)  J 


=  arctan 


sin  (f>G 


cos  <j)G  +  |G| 


(3.12) 


For  ^cl=— 90°  (defines  closed-loop  bandwidth),  the  following  two  equations  must  be  met. 


sin  (j)G  <  0 


(3.13) 


cos  <f>G  +  \G\  =  0 


(3.14) 


Equation  (3.14)  can  be  used  to  solve  for  K  from  Eq.  (3.8). 


|G|  =  \K\  ■  \YVYC\  =  -  cos  <j)G 


\K\ 


-  COS  (f)Q 

W 


In  order  for  \G\  >0, 


COS  (j)G  <  0. 


(3.15) 

(3.16) 

(3.17) 


The  only  place  where  both  Eq.  (3.13)  and  Eq.  (3.17)  are  satisfied  is  when 


-180°  <  <j)G  <  -90°. 


(3.18) 
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If  all  the  parameters  that  make  up  YPYC  are  known  except  for  TPl  and  TP2,  | YPYC\  and  <\>g 
become  functions  of  TPl  and  TP2. 

\YpYc\  =  f1(Tin,Tpa)  and  cf>G  =  f2(TPl,TP2).  (3.19) 

Therefore,  K  is  a  function  of  TPl  and  TP2.  Equation  (3.16)  is  the  equation  that  reduces 
the  two-dimensional  problem  from  Chapter  2  down  to  a  one-dimensional  search  over  TPl .  If 
‘optimum’  lead  compensation  is  used,  TP2  is  known  (TP2  =  0).  If  ‘optimum’  lag  compensation 
is  known,  TP2  is  a  function  of  only  TPl  (assuming  uBw  is  fixed).  This  makes  K  a  function 
of  only  TPl .  The  original  three-dimensional  problem  is  finally  reduced  to  a  one-dimensional 
search  over  TPl .  Now  the  approach  can  be  outlined. 

3.5  Step-by-Step  Procedure 

This  section  only  applies  to  mapping  the  u>spTg2  ,  (sp,  rg  criteria  into  the  Neal-Smith  and 
Pilot-in-the-Loop  criteria.  Figure  3.2  provides  a  flow  chart  of  a  major  part  of  the  process.  The 
steps  before  the  flow  chart  begins  involve  fixing  some  of  the  parameters  of  Yc  in  Eq.  (3.7). 

Step  1:  Select  a  flight  phase.  This  determines  ojbw-  MIL-STD-1797A  lists  two 
different  values  of  u>Bw  for  Category  C  (4:237).  For  landing,  toBW  =  2.5  (rad/sec).  For  any 
other  Category  C  maneuver,  uBw  =  1.5  (rad/sec). 

Step  2 :  Fix  t6  .  Depending  on  whether  a  level  1  region  or  a  level  2  region  is  to  be  mapped 
determines  the  maximum  allowable  rg.  Table  2.1  lists  these  limits.  The  most  conservative 
mapping  would  fix  rg  at  one  of  the  maximum  values  (e.g.  t0  =  0.10  seconds  for  level  1). 

Step  3:  Fix  Tg2.  Choosing  an  aircraft  fixes  Tg2.  For  example,  the  F-16  VISTA  has 
Tg2  =  1.96. 

Step  4:  Check  if  Tg2  is  within  correct  limits.  Table  2.2  outlines  these  limits. 

Step  5:  Determine  the  minimum  value  of  wspTg2.  Table  2.2  also  lists  minimum  values 
for  losp.  Multiplying  one  of  these  values  by  Tg2  determines  the  minimum  value  of  tospTg2.  For 
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that  achieves 
-90°  C-L  phase 

that  achieves 
-90°  C-L  phase 

rc  <T  droop  >  -3  dB^>-^- 


calculate  droop 


droop  =  -3  dB 


calculate  droop 


droop  =  -3  dB 


Figure  3.2  Process  Flow  Chart 


3-8 


example,  the  minimum  value  of  cospTg2  for  a  level  1  region  for  the  F-16  VISTA  would  be  1.7. 
Note  that  this  is  above  the  1.3  boundary  shown  in  Fig.  2.1. 

Step  6:  Determine  range  for  usp  and  (sp  that  defines  a  region  to  be  mapped.  If  the  level 
1  region  is  to  be  mapped,  Eq.  (3.20)  and  Eq.  (3.21)  would  be  the  range. 


10 

0.87  <  c osp<  — 

1$2 

(3.20) 

0.35  <  (sp  <  1.30 

(3.21) 

It  should  be  noted  there  is  no  upper  bound  on  cospTg2.  Since  MIL-STD-1797A  only  plots  the 
criteria  up  to  ijospTg2  =  10  (4:191),  the  upper  bound  in  Eq.  (3.20)  will  be  used. 

Step  7:  Choose  one  point  that  satisfies  Eq.  (3.20)  and  Eq.  (3.21).  At  this  point,  all  the 
parameters  of  Yc  are  known. 

Step  8:  Calculate  the  phase  of  Yc  and  determine  whether  lead  or  lag  compensation  is 
required.  This  is  where  the  flow  chart  in  Fig.  3.2  begins.  If  LYC  <  -180°  or  LYC  is  in  the  target 
zone  but  droop  <  -3  dB,  lead  compensation  is  needed.  ‘Optimal’  lead  compensation  will  be 
used,  so  Eq.  (2.18)  applies.  If  LYC  >  —125.3°  or  LYC  is  in  the  target  zone  but  droop  >  -3  dB, 
lag  compensation  is  needed.  ‘Optimal’  lag  compensation  will  be  used,  so  Eq.  (2.20)  applies. 

There  are  limits  on  the  amount  of  phase  the  pilot  model  can  add  or  subtract.  Figures  2.1 1 
and  2.16  asymptotically  approach  ±90°  respectively.  Therefore,  the  maximum  amount  of 
phase  the  pilot  model  can  add  or  subtract  is  90°.  If  closed-loop  bandwidth  and  droop  are  to 
be  met  simultaneously,  Eq.  (3.22)  must  be  satisfied.  The  computer  code  is  designed  to  check 
for  this  and  flag  those  points  that  do  not  satisfy  Eq.  (3.22). 

-  270°  <  LYC  <  -35.3°  (3.22) 

Step  9:  Find  ‘optimal’  values  for  K,  TPl,  and  TP2  by  doing  a  one-dimensional  search 
over  TPl.  Here,  ‘optimal’  means  a  combination  of  K,  TPl,  and  TP2  that  makes  the  pilot- 
aircraft  system  meet  the  closed-loop  bandwidth  and  droop  requirements  simultaneously.  The 
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command  ‘constr’  in  MATLAB™  solves  constrained  optimization  problems  using  SQP. 
Tolerances  are  needed  to  determine  if  convergence  is  achieved.  The  default  tolerances  in  the 
Optimization  Toolbox  (8)  are  used  and  shown  in  Table  3.3. 


Table  3 .3  Tolerances  for  Convergence 


Parameter 

Tolerance 

Objective  function  (f) 

10“4 

Constraint  violation  ( g ) 

10~6 

Design  vector  (x) 

io-4 

Droop  is  determined  by  doing  a  frequency  search  for  the  minimum  value  of  \T  \  over  lo  <  ojbw- 
The  ‘min’  command  in  MATLAB™  accomplishes  this.  The  design  vector  becomes  a  scalar 
( TP1 ).  The  cost  function  becomes  the  square  of  the  difference  between  the  closed-loop  droop 
and  -3  dB.  The  specific  problem  statement  for  optimizing  with  lead  compensation  is  given  by 
Table  3.4.  The  specific  problem  statement  for  optimizing  with  lag  compensation  is  given  by 
Table  3.5. 

Side  constraints  are  placed  on  the  pilot  parameters  (TPl  >  0  and  K  >  0).  Negative 
values  would  lead  to  non-minimum  phase  systems.  The  next  section  outlines  the  approach  for 
identifying  regions  of  agreement  and  regions  of  conflict  after  the  mapping  of  u>spTg2,(Sp,Te 
is  complete. 

Table  3.4  Lead  Problem  Statement 


Minimize:  ( DROOP  +  3)2  (dB)2 

Subject  to:  DROOP  >  —3  dB 

-180°  <<j>G<  -125.3° 
K  >  0 
TP1  >0 
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Table  3.5  Lag  Problem  Statement 


Minimize:  ( DROOP  +  3)2  (dB)2 

Subject  to:  DROOP  >  —3  dB 

-180°  <  4>g  <  -125.3° 
K  >  0 
TPi  >0 
TP2  >  TP1 


3.6  Region  Identification 

Both  the  Pilot-in-the-Loop  and  Neal-Smith  criteria  involve  finding  values  for  the  pilot 
model  with  respect  to  a  given  aircraft  model.  A  point  from  either  criteria  represents  the 
pilot-aircraft  system  in  the  form  of  one  transfer  function.  The  individual  pilot  and  aircraft 
parameters  are  not  easily  distinguished.  Rather  than  trying  to  accomplish  this  for  a  grid  of 
points,  bookkeeping  of  the  original  mapping  from  the  last  section  will  be  used.  Entire  regions 
from  either  criteria  will  never  be  mapped.  Only  those  parts  corresponding  to  the  region 
from  uspTo2,(Sp,Te  will  be  evaluated.  The  identification  of  conflict  and  agreement  regions 
involves  pass-fail  logic.  For  example,  when  a  level  1  point  from  lo3PT$2,(sP,  tq  is  mapped 
into  Neal-Smith,  does  it  fall  within  the  Neal-Smith  level  1  boundaries?  If  the  answer  is  yes,  a 
certain  flag  is  attached  to  that  point.  If  the  answer  is  no,  a  different  flag  is  used.  All  the  points 
with  the  ‘yes’  flag  can  be  mapped  back  to  the  uspTg2,(sp,  rg  criteria  showing  the  region  of 
agreement.  The  region  of  conflict  is  identified  in  the  same  fashion.  Before  beginning  Chapter 
IV,  a  summary  of  Chapter  III  is  given. 

3.7  Summary 

The  original  three-dimensional  problem  of  determining  K,  TPl,  and  TP2  was  reduced 
to  a  one-dimensional  search  over  TPl.  With  all  of  the  parameters  fixed  except  TP1  and 
TP2,  the  open-loop  transfer  function  G  is  a  function  of  only  TPl  and  TV2.  If  TP1  is  known, 
TP2  =  ■  TP1 )  when  ‘optimum’  lag  compensation  is  used  and  TP2  =  0  when  ‘optimum’ 
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lead  compensation  is  used.  Therefore  G  =  j\{TPl).  A  formula  for  \I<\  as  a  function  of  the 
magnitude  and  phase  of  G  results  in  K  =  f2 (G)  =>•  I(  =  f2(TPl).  Hence,  if  TPl  is  known,  so 
is  I<  and  TP2.  The  closed-loop  bandwidth  and  droop  requirements  are  turned  into  constraints. 
If  LYC  is  to  left  of  the  target  zone,  lead  compensation  is  needed.  If  lYc  is  to  right  of  the  target 
zone,  lag  compensation  is  needed.  If  lYc  is  in  the  target  zone,  lead  compensation  is  needed  if 
the  droop>-3  dB  and  lag  compensation  is  needed  if  the  droop<-3  dB.  Both  the  lead  and  lag 
problem  fit  the  general  form  for  SQP  to  solve.  Once  a  solution  for  the  pilot  model  is  found, 
the  mapping  of  a  region  from  c ospTe2 ,  (sp,  tq  into  either  Neal-Smith  or  Pilot-in-the-Loop  can 
be  accomplished.  From  there,  bookkeeping  of  the  original  mapping  produces  the  ‘inverse’ 
mapping. 
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IV.  Results 


In  order  to  perform  a  mapping,  the  parameters  lobw,  t$,  and  rl\  need  to  be  fixed.  There 
is  an  infinite  number  of  combinations.  Twelve  combinations  are  given  in  this  thesis.  The 
intent  is  to  provide  a  sample  that  will  illustrate  trends  and  highlight  problems. 

Since  lobw,  n,  and  Tg2  are  fixed,  trends  can  be  examined  when  different  values  are 
used.  For  example,  ooBw  =  2.5  (rad/sec)  for  landing  and  lobw  =  1-5  (rad/sec)  for  any 
other  Category  C  maneuver.  Holding  to  and  Tg,2  constant  while  lobw  varies  enables  the  trend 
of  increasing  lobw  to  be  examined.  Holding  lobw  and  Tg2  constant  enables  the  trend  of 
increasing  tq  to  be  examined,  tq  will  be  placed  at  its  maximum  value  depending  on  which 
handling  qualities  level  is  being  mapped  (see  Table  2.1).  t6  will  also  be  placed  at  the  level 
1  maximum  value  for  level  2  mappings.  Finally,  holding  lobw  and  Tg  constant  for  the  two 
different  aircraft  enables  the  trend  of  increasing  Tg2  to  be  examined. 

In  total,  six  cases  will  be  examined  for  each  aircraft.  Table  4.1  defines  each  case.  The 
‘Level’  column  refers  to  the  region  from  the  LospTg2 ,  (sp ,  Tg  criteria  that  is  mapped.  The  region 
is  formed  by  fixing  u>bw,  to,  and  Tg2  to  their  respective  values  while  varying  cosp  and  £ sp . 
The  boundaries  defined  in  Fig.  2.1  are  used.  However,  the  minimum  boundary  for  uspTg2 
is  different  for  the  F-16  VISTA.  The  minimum  values  of  uosp  in  Table  2.2  come  into  effect. 
The  minimum  boundary  for  level  1  becomes  1.7  instead  of  1.3  as  in  Fig.  2.1.  For  level  2, 
the  minimum  boundary  becomes  1.18  instead  of  0.75.  The  results  of  a  sample  point  will 
be  presented  followed  by  the  results  of  a  sample  case.  The  results  of  all  12  cases  are  in 
Appendices  (A-L). 

4.1  Sample  Point 

This  section  will  show  the  information  required  to  map  one  point  in  case  2  from  the 
iospTg2 ,  £sp,  Tg  criteria  into  both  the  Neal-Smith  and  Pilot-in-the-Loop  criteria.  The  point  will 
have  level  1  handling  qualities  for  the  LospTg2,(sp,  Tg  and  the  Pilot-in-the-Loop  criteria,  but 
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Table  4.1  Test  Cases 


Case 

Level 

UJBW 

TO 

To 2 

1 

1 

1.50 

0.10 

1.67 

2 

1 

2.50 

0.10 

1.67 

3 

2 

1.50 

0.10 

1.67 

4 

2 

2.50 

0.10 

1.67 

5 

2 

1.50 

0.20 

1.67 

6 

2 

7 

1.96 

8 

1 

mi 

9 

mm 

10 

mm 

1.96 

2 

1.50 

0.20 

1.96 

12 

2 

2.50 

0.20 

1.96 

will  have  level  2  handling  qualities  for  the  Neal-Smith  criteria.  The  required  parameters  that 
define  the  point  are  shown  in  Table  4.2. 


Table  4.2  Sample  Point 


Case 

Level 

UBW 

re 

To  2 

^ sp 

Csp 

2 

1 

2.50 

0.10 

1.67 

1.25 

0.83 

Equation  4.1  becomes  the  aircraft  transfer  function. 

(1.67  -5  +  1)  e~010s 
c  “  s{s2  +  2.085  +  1.56) 


(4.1) 


Figure  4.1  illustrates  that  the  sample  point  is  contained  by  the  level  1  boundary  of  the 
u!spTg2,(sp,  tq  criteria.  The  next  step  involves  determining  whether  lead  or  lag  compen¬ 
sation  is  needed.  On  Fig.  4.2,  (o)  marks  the  open-loop  magnitude  and  phase  of  Yc  where 
to  —  lujjw-  Since  LYC  <  —180°  at  u;  =  loBw,  lead  compensation  is  needed.  Figure  4.3 
shows  the  result  of  K  ■  YPYC  after  ‘optimum’  lead  compensation  has  been  applied.  Both  the 
closed-loop  bandwidth  requirement  and  closed-loop  phase  requirement  are  met  simultane¬ 
ously.  From  Fig.  4.4,  the  closed-loop  droop  satisfies  the  MIL-STD-1797A  requirements,  and 
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Figure  4. 1  Sample  Point  on  the  u>spTe2 ,  (sp ,  tq  Criteria 


the  closed-loop  resonance  is  well  below  the  3  dB  line.  Figure  4.5  shows  the  closed-loop  phase 
equals  -90  degrees  when  to  =  ubw-  Hence,  the  point  maps  to  level  1  for  the  Pilot-in-the-Loop 
criteria.  Figure  4.6,  however,  shows  the  pilot’s  phase  is  high  enough  to  map  the  point  to  the 
level  2  region  on  the  Neal-Smith  criteria.  This  sample  point  illustrates  a  conflicting  result 
between  the  Neal-Smith  and  Pilot-in-the-Loop  criteria.  A  sample  case  is  mapped  in  the  next 
section.  Conflicts  among  all  three  criteria  will  be  illustrated. 

The  computer  code  allows  Fig.  4.1  -  Fig.  4.6  to  be  viewed  simultaneously  on  a  computer 
monitor.  When  mapping  a  grid  of  points,  the  monitor  displays  all  six  plots  for  each  point  in 
a  ‘picture  show’  fashion.  In  other  words,  the  monitor  displays  the  results  frame  by  frame  at 
a  user  specified  time  increment.  This  tool  enables  the  user  to  see  trends  as  they  occur  and 
ensure  the  closed-loop  bandwidth  and  droop  requirements  are  being  met  for  each  point. 
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Open-Loop  Mag  dB  Open-Loop  Mag  dB 


Figure  4.2  Nichols  Chart  of  Yc 


Figure  4.3  Nichols  Chart  of  YCYP  after  Convergence 
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Closed-Loop  Phase  (de 


Figure  4.6  Sample  Point  Mapped  to  Neal-Smith  Boundaries 


4.2  Sample  Case 

This  section  will  show  all  the  plots  involved  in  mapping  1  of  the  12  cases.  Case  2  from 
Table  4.1  will  be  mapped  because  it  shows  that  infeasible  regions  are  possible.  Here,  infeasible 
means  the  closed-loop  bandwidth  and  droop  requirements  cannot  be  met  simultaneously. 
When  this  happens,  the  ‘constr’  command  in  MATLAB™  produces  a  ‘Warning:  No  feasible 
solution  found’  message.  Figure  4.7  shows  the  Nichol’s  plot  of  Yc  for  one  of  the  points  in 
the  infeasible  region  shown  in  Fig.  4.8.  The  point  where  lo  =  u>bw  (°)  is  to  the  left  of  the 
target  zone,  and  the  magnitude  is  0  dB.  Also,  the  droop  is  below  the  -3  dB  closed-loop  droop 
line.  The  —90°  closed-loop  phase  line  does  not  extend  above  the  0  dB  open-loop  magnitude 
line.  Lead  compensation  would  cause  the  point  to  move  up  and  to  the  right.  The  open-loop 
magnitude  of  (o)  would  move  above  the  0  dB  line.  The  gain  could  be  adjusted  to  lower 
the  curve  to  meet  the  closed-loop  bandwidth  constraint.  This  would  cause  the  droop  to  fall 
back  below  the  -3  dB  closed-loop  droop  line.  The  result  is  both  constraints  cannot  be  met 
simultaneously.  Hence,  the  point  cannot  be  mapped  into  the  Neal-Smith  criteria  and  will 
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Figure  4.7  Infeasible  Point  from  the  LospTe2 ,  (sp,  tb  Criteria 


be  considered  infeasible.  Infeasible  regions  occur  because  the  assumed  pilot  model  cannot 
provide  the  required  compensation  needed  to  meet  both  constraints.  These  regions  occur  for 
both  lead  and  lag  compensation.  The  results  show  that  infeasible  regions  are  only  found  in 
regions  corresponding  to  Neal-Smith  and  Pilot-in-the-Loop  level  3.  The  next  few  paragraphs 
describe  the  plots  involved  in  mapping  Case  2.  Six  plots  will  be  made  for  each  case.  Again, 
the  plots  for  all  12  cases  are  shown  in  Appendices  (A-L). 

4.2.1  Plot  1.  This  plot  shows  the  region  from  the  wspTe2 ,  (sp,  tb  criteria  that  will 
be  mapped.  The  region  is  divided  into  areas  requiring  lead  compensation  and  areas  requiring 
lag  compensation.  The  infeasible  area  is  also  shown.  Figure  4.8  shows  that  the  dividing  line 
between  lead  and  lag  compensation  is  fairly  smooth.  If  the  logic  from  Chapter  III  determines 
that  lead  compensation  is  needed,  that  point  is  flagged  ‘lead’.  Otherwise,  it  is  flagged  ‘lag’. 
The  original  grid  can  then  be  divided  into  a  lead  area  and  a  lag  area.  Corners  on  the  grid, 
as  well  as  other  dividing  points,  are  marked  with  letters  to  show  where  those  points  end  up 
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Figure  4.8  (Plot  1)  Lead,  Lag,  and  Infeasible  Regions  on  uspT62 ,  (sp,  t6 

after  the  mapping.  Level  1  mappings  use  an  evenly  spaced  30x50  grid  on  a  log-log  scale. 
The  approach  in  Chapter  III  is  called  1500  times.  Data  required  to  produce  the  information 
in  Fig.  4.1  -  Fig.  4.6  is  calculated  and  stored  for  all  1500  points.  A  mapping  of  this  size  takes 
about  three  hours  to  run  on  a  Sun  SPARCstation  20.  Level  2  mappings  use  an  evenly  spaced 
50x60  grid  on  a  log-log  scale.  The  computer  time  doubles  by  increasing  the  number  of  points 
to  3000. 

4.2.2  Plot  2.  This  plot  shows  the  region  from  plot  1  mapped  into  the  Neal-Smith 
criteria.  Figure  4.9  shows  that  conflicts  indeed  exist  between  the  ojspTe2 ,  (sp,  rg  and  Neal- 
Smith  Criteria.  The  shaded  regions  in  Fig.  4.9  were  made  by  connecting  the  outer  points 
that  form  the  boundary  for  a  given  area.  The  marked  points  from  Fig.  4.8  are  also  shown  on 
Fig.  4.9.  The  letters  G,  H,  and  I  are  on  the  border  of  the  infeasible  region  shown  in  Fig.  4.8. 
These  letters  map  well  into  the  Neal-Smith  level  3  region  on  Fig.  4.9. 


4-8 


60 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


-80  -60  -40  -20  0  20  40  60  80  100 


Pilot  Phase  Angle  at  Bandwidth  (deg) 


Figure  4.9  (Plot  2)  iospTg2 ,  (sp,  t6  Mapped  into  Neal-Smith 

4.2.3  Plot  3  and  Plot  4.  Plot  3  is  the  same  as  plot  2  except  regions  corresponding 
to  Neal-Smith  levels  1-3  are  highlighted.  Figure  4.10  shows  the  three  different  handling 
qualities  regions  according  to  Neal-Smith  standards.  Plot  4  shows  these  highlighted  regions 
mapped  back  into  the  cospTg2,(sp,  t9  criteria.  Like  Fig.  4.9,  Fig.  4.11  shows  both  regions  of 
agreement  and  regions  of  conflict.  The  region  of  agreement  marked  by  J,  K,  and  L  is  fairly 
small  compared  to  original  level  1  region  from  the  uspTg2 ,  (sp,  t8  criteria. 

4.2.4  Plot  5  and  Plot  6.  Plot  5  is  the  same  as  plot  2  except  regions  corresponding 
to  Pilot-in-the-Loop  levels  1-3  are  highlighted.  Figure  4.12  shows  the  three  different  handling 
qualities  regions  according  to  Pilot-in-the-Loop  standards.  Plot  6  shows  these  highlighted 
regions  mapped  back  into  the  uspTg2,(sp,  t0  criteria.  Not  surprisingly,  the  level  1  region  of 
agreement  between  the  uspTg2,(sp,  Tg  and  Pilot-in-the-Loop  criteria  is  larger  than  the  level 
1  region  of  agreement  between  the  LospTg2,(sP ,  Tg  and  Neal-Smith  criteria.  The  Pilot-in-the- 
Loop  criteria  does  not  take  pilot  phase  angle  into  consideration.  Hence,  its  level  1  boundary 
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Pilot  Phase  Angle  at  Bandwidth  (deg) 


Figure  4. 1 1  (Plot  4)  Neal-Smith  Regions  Mapped  into  a >spTe2 ,  ( 
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Figure  4. 12  (Plot  5)  Pilot-in-the-Loop  Regions  (PL) 

shows  up  a  straight  line  at  3  dB  on  Fig.  4.10.  Compared  to  the  Neal-Smith  level  1  boundary  in 
Fig.  4.9,  the  Pilot-in-the-Loop  will  always  cover  an  area  greater  than  or  equal  to  the  Neal-Smith 
criteria.  This  is  evident  when  comparing  Fig.  4.11  to  Fig.  4.13. 

The  level  2  region  is  not  as  straightforward.  The  Pilot-in-the-Loop  level  2  border  is 
a  straight  line  at  9  dB  on  Fig.  4.10.  The  Neal-Smith  level  2  border  is  greater  than  9  dB  for 
pilot  phase  angles  less  than  about  40°  and  less  than  9  dB  elsewhere.  The  level  1  borders  for 
both  criteria  have  already  been  shown  to  be  different.  The  level  1  borders  determine  the  lower 
limit  of  the  level  2  regions.  Since  both  the  upper  and  lower  limits  are  different,  no  general 
conclusions  can  be  made.  For  the  sample  case,  the  region  corresponding  to  Pilot-in-the-Loop 
level  3  was  greater  than  Neal-Smith  level  3.  This  again  can  be  seen  by  comparing  Fig.  4.1 1  to 
Fig.  4.13. 
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Figure  4.13  (Plot  6)  Pilot-in-the-Loop  Mapped  into  the  uspTg2 ,  (sp,  Te 
4.3  Regions  Excluded  by  (u>sp )min 

As  explained  earlier,  the  minimum  boundary  for  iospTg2  is  moved  up  for  the  F-16  VISTA. 
This  excludes  parts  of  the  general  level  1  and  level  2  regions  on  the  u>spTg2,(sp,  t0  criteria. 
These  regions  were  mapped  separately  for  Cases  7-12.  Plots  4  and  6  were  also  generated  for 
the  excluded  regions  and  appear  in  Appendices  (G-L).  They  are  designated  with  an  ‘x\  Since 
the  excluded  regions  are  fairly  small,  only  one  plot  is  generated  to  show  the  uspTg2,(sP,  Tg 
criteria  mapping.  Both  Neal-Smith  and  Pilot-in-the-Loop  boundaries  are  shown  on  the  same 
plot.  Individual  points  are  used  instead  of  the  shaded  regions.  Figure  4.14  is  a  sample  of  this 
mapping.  As  seen  in  Appendices(G-L),  the  excluded  regions  map  to  areas  in  the  lower  right 
corner  of  the  Neal-Smith  criteria.  This  generates  a  conflict  for  level  1  mappings. 

The  Pilot-in-the-Loop  criteria  stayed  one  level  above  the  Neal-Smith  criteria.  If  the 
excluded  region  mapped  to  Neal-Smith  level  2,  it  mapped  to  Pilot-in-the-Loop  level  1.  Again, 
this  is  not  surprising  since  the  Pilot-in-the-Loop  criteria  does  not  account  for  pilot  phase.  The 
only  general  statement  that  can  be  made  is  the  same  as  above.  The  Neal-Smith  criteria  will 
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Figure  4.14  (Case  7x)  Mapping  of  Excluded  Region  from  LospTg2 ,  (sp,  re 


never  predict  level  1  handling  qualities  more  times  than  the  Pilot-in-the-Loop  criteria.  That 
is,  the  Pilot-in-the-Loop  criteria  is  more  liberal  for  level  1  handling  qualities. 


4.4  Trends 

For  a  better  understanding  of  trends,  all  the  plots  in  Appendices  (A-L)  need  to  be 
examined.  This  section  will  present  a  simplified  analysis  of  trends.  Rather  than  plotting  all 
12  cases  on  a  plot,  an  example  comparing  two  cases  is  shown  for  each  trend.  The  examples 
represent  general  trends.  Trends  shown  are  for  the  u>spTg2 ,  (sp ,  t0  and  Neal-Smith  mappings. 

4.4.1  Effect  of  Changing  c oBw-  The  figures  referenced  in  this  paragraph  apply  to  the 
Lear  jet.  Increasing  ojbw  from  1.5  to  2.5  (rad/sec)  causes  the  line  dividing  lead  and  lag  regions 
on  the  UspTe2,(sp,Te  criteria  to  move  up  as  shown  in  Fig.  4.15.  The  region  corresponding 
to  Neal-Smith  and  uspTg2,(sp,  rg  level  1  moves  up  and  to  the  right  as  shown  in  Fig.  4.16. 
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Finally,  the  region  corresponding  to  Neal-Smith  level  2  moves  up  and  to  the  right  as  shown  in 
Fig.  4.17. 

4.4.2  Effect  of  Changing  Tg2 .  The  figures  referenced  in  this  paragraph  correspond 
to  lobw  =  1-5  (rad/sec).  Increasing  Tg2  from  1.67  (Lear  jet)  to  1.96  (F-16  VISTA)  causes 
the  line  dividing  lead  and  lag  regions  on  the  LospTg2,(sp,  t6  criteria  to  move  up  as  shown  in 
Fig.  4.18.  The  region  corresponding  to  Neal-Smith  and  eospTg2 ,  (sp,  r6  level  1  moves  slightly 
up  and  to  the  right  as  shown  in  Fig.  4.19.  Finally,  the  region  corresponding  to  Neal-Smith 
level  2  moves  up  and  to  the  right  as  shown  in  Fig.  4.20. 

4.4.3  Effect  of  Changing  t0.  The  figures  referenced  in  this  paragraph  correspond 
to  lobw  =  1.5  (rad/sec).  Increasing  tq  from  0.1  to  0.2  seconds  causes  the  line  dividing  lead 
and  lag  regions  on  the  u)spTg2,(sp,  Tg  criteria  to  move  up  as  shown  in  Fig.  4.21.  The  region 
corresponding  to  Neal-Smith  and  iospTg2 ,  (sp,  tb  level  1  moves  up  and  to  the  right  as  shown  in 
Fig.  4.22.  Finally,  the  region  corresponding  to  Neal-Smith  level  2  moves  up  and  to  the  right 
as  shown  in  Fig.  4.23. 

4.5  Summary 

The  information  required  to  map  12  cases  was  outlined  in  Table  4.1.  The  results  of 
a  sample  point  and  sample  case  were  discussed.  Parametric  studies  were  performed  and 
analyzed  for  trends.  Increasing  luBw,  Tg2,  or  re  resulted  in  smaller  regions  of  agreement. 
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Figure  4. 15  Effect  of  Increasing  u>bw  on  Lead  and  Lag 


Figure  4.16  Effect  of  Increasing  uBW  on  Level  1  Region  of  Agreement 


10 


Figure  4.19  Effect  of  Increasing  Tg2  on  Level  1  Region  of  Agre 


Figure  4.20  Effect  of  Increasing  Tg2  on  NS  Level  2  Regior 
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Figure  4.21  Effect  of  Increasing  tq  on  Lead  and  Lag 


Figure  4.22  Effect  of  Increasing  re  on  Level  1  Region  of  Agreement 


'  10° 


Figure  4.23  Effect  of  Increasing  tq  on  NS  Level  2  Region 
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V.  Conclusions  and  Recommendations 


5.1  Conclusions 

All  the  objectives  for  this  thesis  were  accomplished.  The  development  of  a  computerized 
tool  that  performs  a  mapping  of  the  c ospTo2 ,  (sp,  t8  criteria  into  the  Neal-Smith  criteria  has  been 
presented.  The  problem  of  calculating  the  parameters  of  the  Neal-Smith  pilot  model  was  posed 
as  a  constrained  optimization  problem.  Initially,  the  problem  involved  the  three  unknowns  K, 
Tp  15  and  TP2.  Deriving  equations  for  ‘optimal’  lead  and  ‘optimal’  lag  compensation  left  TP2 
in  terms  of  TPl.  Deriving  an  equation  for  the  closed-loop  phase  and  setting  it  equal  to  —90° 
left  K  in  terms  of  TVl .  Finding  the  solution  that  met  both  the  closed-loop  bandwidth  and 
droop  requirements  came  down  to  a  one-dimensional  search  over  TPl .  Sequential  Quadratic 
Programming  was  needed  to  solve  for  TPl  because  no  closed  form  solutions  exist  for  the 
complicated  functions  of  frequency.  Bookkeeping  of  the  original  mapping  produced  the 
information  required  to  do  an  inverse  mapping. 

Six  cases  for  two  Tg2  values  were  mapped.  The  u>spTg2,(sp,  t8  criteria  was  mapped 
into  both  the  Neal-Smith  and  Pilot-in-the-Loop  criteria.  The  regions  corresponding  to  both 
Neal-Smith  and  Pilot-in-the-Loop  levels  1-3  were  mapped  back  into  uispTg2,(sp,  re  criteria. 
Regions  of  agreement  and  regions  of  conflict  could  be  easily  seen  in  the  six  plots  that  were 
generated  for  each  case.  The  trends  of  varying  the  three  fixed  parameters  were  identified. 

5.2  Recommendations  for  Future  Research 

The  results  produced  by  this  thesis  lay  the  foundation  for  a  flight  test  program.  Selected 
points  from  any  region  of  conflict  could  be  flown  to  determine  the  actual  handling  qualities. 
The  results  could  then  be  analyzed  to  determine  if  either  the  a )spTg2 ,  (sp,  t6  or  Neal-Smith 
criteria  needs  modification.  The  information  learned  from  this  type  of  analysis  is  the  next  step 
towards  deriving  one  clear-cut,  generally  applicable  set  of  requirements. 


5-1 


Appendix  A .  Case  1 
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Figure  A.l  (Case  1)  Lead,  Lag,  and  Infeasable  Regions  on  &SpTg 2 ,  Cspi 


Figure  A.2  (Case  1 )  cosp  Tg2 ,  (sp ,  t0  Mapped  into  Neal-Smith 
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60 


Figure  A.3  (Case  1)  Neal-Smith  Regions  (NS) 


Figure  A.4  (Case  1)  Neal-Smith  Regions  Mapped  into  ojspTg2 ,  (sp ,  ts 
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Appendix  B.  Case  2 
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Figure  B.l  (Case  2)  Lead,  Lag,  and  Infeasible  Regions  on  cospTe2 ,  (sp 
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Figure  B.2  (Case  2)  uspTq2  ,  (sp ,  t0  Mapped  into  Neal-Smith 
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Figure  B.3  (Case  2)  Neal-Smith  Regions  (NS) 


Figure  B.4  (Case  2)  Neal-Smith  Regions  Mapped  into  wspTg2 ,  (sp,  t6 
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Appendix  C.  Case  3 
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Figure  C.2  (Case  3)  uspTg2 ,  (sp,  re  Mapped  into  Neal-Smith 
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Figure  C.4  (Case  3)  Neal-Smith  Regions  Mapped  into  uspTe2 ,  (sp,  t6 


Figure  C.5  (Case  3)  Pilot-in-the-Loop  Regions  (PL) 


Figure  C.6  (Case  3)  Pilot-in-the-Loop  Mapped  into  the  ujspTe2 ,  (sp ,  re 


Appendix  D.  Case  4 
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Figure  D.2  (Case  4)  uspTe2 ,  (sp,  t6  Mapped  into  Neal-Smith  • 
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Pilot  Phase  Angle  at  Bandwidth  (deg) 


Figure  D.4  (Case  4)  Neal-Smith  Regions  Mapped  into  uapTg2 ,  (sp,  re 


Appendix  E.  Case  5 
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Figure  E.2  (Case  5)  u}spTg2 ,  (sp ,  r@  Mapped  into  Neal-Smith 


Pilot  Phase  at  Bandwidth  (deg) 


Figure  E.4  (Case  5)  Neal-Smith  Regions  Mapped  into  tospTg2 ,  (sp,  rg 


Figure  E.5  (Case  5)  Pilot-in-the-Loop  Regions  (PL) 


Appendix  F.  Case  6 
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Figure  F.2 
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(Case  6)  uspTg2,(sp, t6  Mapped  into  Neal-Smith 
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Figure  F.3  (Case  6)  Neal-Smith  Regions  (NS) 
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Figure  F.4  (Case  6)  Neal-Smith  Regions  Mapped  into  u}spTg2 ,  (sp,  re 
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Pilot  Phase  at  Bandwidth  (deg) 


Figure  F.5  (Case  6)  Pilot-in-the-Loop  Regions  (PL) 
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Figure  F.6  (Case  6)  Pilot-in-the-Loop  Mapped  into  the  iospTg2 ,  (sp,  tb 


Appendix  G.  Case  7 
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Figure  G.l  (Case  7)  Lead,  Lag,  and  Infeasible  Regions  on  u>spTe2 ,  (sp,  re 
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Figure  G.2  (Case  7)  uspTo2 ,  (sp,  re  Mapped  into  Neal-Smith 
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Figure  G.4  (Case  7)  Neal-Smith  Regions  Mapped  into  wspTo2 ,  £ Sp  j  re 


Figure  G.5  (Case  7)  Pilot-in-the-Loop  Regions  (PL) 


Figure  G.6  (Case  7)  Pilot-in-the-Loop  Mapped  into  the  uspTe2 ,  (sp ,  re 
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Appendix  H.  Case  8 
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Figure  H.2  (Case  8)  usp  Tg2 ,  (sp ,  ts  Mapped  into  Neal-Smith 
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Figure  H.5  (Case  8)  Pilot-in-the-Loop  Regions  (PL) 
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Figure  H.7  (Case  8x)  Mapping  of  Excluded  Region  from  cospTg2 ,  (sp ,  tb 


Figure  H.8  (Case  8x)  Neal-Smith  Regions  Mapped  into  u:spTg2 ,  (sp ,  t0 
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Appendix  I.  Case  9 
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Figure  1.2  (Case  9)  uspTo2 ,  (sp,  re  Mapped  into  Neal-Smith 
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Figure  1.3  (Case  9)  Neal-Smith  Regions  (NS) 


Figure  1.4  (Case  9)  Neal-Smith  Regions  Mapped  into  ojspTe2 ,  (sp,  t6 
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Figure  1.5  (Case  9)  Pilot-in-the-Loop  Regions  (PL) 


Figure  1.6  (Case  9)  Pilot-in-the-Loop  Mapped  into  the  ojspTg2 ,  (spi  t6 
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Appendix  J.  Case  10 
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Figure  J.  1  (Case  10)  Lead,  Lag,  and  Infeasible  Regions  on  u>spT$2 ,  (sp ,  re 
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Figure  J.2  (Case  10)  t ospTe2 ,  (sp,  re  Mapped  into  Neal-Smith 
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Figure  J.3  (Case  10)  Neal-Smith  Regions  (NS) 


Figure  J.4  (Case  10)  Neal-Smith  Regions  Mapped  into  u)spTg2 ,  (sp ,  t0 
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Figure  J.5  (Case  10)  Pilot-in-the-Loop  Regions  (PL) 


Figure  J.6  (Case  10)  Pilot-in-the-Loop  Mapped  into  the  uspT62 ,  (sp ,  t0 
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Figure  J.8  (Case  lOx)  Neal-Smith  Regions  Mapped  into  u>spTg2 ,  (sp ,  t6 


Appendix  K.  Case  11 
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Figure  K.2  (Case  11)  u)sp  Tg2 ,  (sp ,  rg  Mapped  into  Neal-Smith 
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Figure  K.5  (Case  11)  Pilot-in-the-Loop  Regions  (PL) 
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Figure  K.6  (Case  1 1)  Pilot-in-the-Loop  Mapped  into  the  wspTg2 ,  (sp ,  re 
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Figure  K.7  (Case  1 1  x)  Mapping  of  Excluded  Region  from  usp  T&2 ,  (sp ,  t0 


Figure  K.8  (Case  1  lx)  Neal-Smith  Regions  Mapped  into  uspTg2 ,  (sp ,  to 


Appendix  L.  Case  12 
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Figure  L.l  (Case  12)  Lead,  Lag,  and  Infeasible  Regions  on  oospT&2  j  Cspi  T 6 
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Figure  L.2  (Case  12)  t ospTe2 ,  (sp,  n  Mapped  into  Neal-Smith 
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Pilot  Phase  at  Bandwidth  (deg) 


Figure  L.3  (Case  12)  Neal-Smith  Regions  (NS) 


Figure  L.4  (Case  12)  Neal-Smith  Regions  Mapped  into  LospTg2 ,  (.sp,  t9 
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Figure  L.5  (Case  12)  Pilot-in-the-Loop  Regions  (PL) 


Figure  L.6  (Case  12)  Pilot-in-the-Loop  Mapped  into  the  ujsvTo2  ,  (sp ,  re 
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Figure  L.8  (Case  12x)  Neal-Smith  Regions  Mapped  into  ojspTg2 ,  (sp,  t6 
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